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Abstract 

In this paper we analyse the vacuum polarization effects associated with a 
massless fermionic field in a higher-dimensional global monopole spacetime in the 
"braneworld" scenario. In this context we admit that the our Universe, the bulk, is 
represented by a flat (n — 1)— dimensional brane having a global monopole in a extra 
transverse three dimensional submanifold. We explicitly calculate the renormalized 
vacuum average of the energy-momentum tensor, (T^(x)}ji en ,, admitting the global 
monopole as being a point-like object. We observe that this quantity depends cru- 
cially on the value of n, and we provide explicit expressions to it for specific values 
attributed to n. 
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1. Introduction 



Recently the braneworld model has attracted renewed interest. By this scenario our world 
is represented by a four dimensional sub-manifold, a three-brane, embedded in a higher 
dimensional spacetime pp. Braneworlds naturally appear in the string/M theory con- 
text and provide a novel setting for discussing phenomenological and cosmological issues 
related to extra dimensions. The models introduced by Randall and Sundrum are partic- 
ularly attractive [21 [3]. The corresponding spacetime contains two (RSI), respectively one 
(RSII), Ricci-flat brane(s) embedded on a five-dimensional Anti-de Sitter (AdS) bulk. It is 
assumed that all matter fields are confined on the branes and only the gravity propagates 
in the five dimensional bulk. The idea that matter is confined to a lower dimensional man- 
ifold is not a new one. The localization of fermions on a domain wall has been discussed 
in [1]. 

The hierarchy problem between the Planck scale and the electroweak one is solved in 
the RSI model, if the distance between the two branes is about 37 times the AdS radius. 
The braneworld model also provide some alternative discussions about one of the most 
important problem in the modern physics: the cosmological constant problem (see, for 
instance, Ref. [5]). In this way, the Casimir energy associated with quantum fields which 
propagates in the bulk obeying specific boundary conditions on the branes may contribute 
to both, the brane and bulk cosmological constant. The Casimir energy associated with 
scalar field on the five-dimensional Randal- Sundrum model are calculated in [Bj. Surface 
Casimir densities and induced cosmological constant on the branes are calculated in [7] 
for a massive scalar quantum field obeying Robin boundary conditions on two parallel 
brane in a general (D + 1)— dimensional anti-de Sitter bulk. 

Although topological defects have been first studied in a four-dimensional space- 
time [8], they have been considered in spacetimes of higher dimensions in the context 
of braneworld. In this scenario the defects live in a d— dimensions submanifold, hav- 
ing their cores on the 3— brane. In this way, domain wall and cosmic string cases have 
been analyzed in [I] and [HUD], respectively, considering d = 1 and 2. Local and global 
monopoles have also been analyzed in [TTJ [12] and [13] to [17], respectively, considering 
d = 3. Specifically in [17J it was shown that if rjo, the energy scale where the gauge 
symmetry of the global system is spontaneously broken, is smaller than the Planck mass, 
the seven-dimensional Einstein equations admit a solution which for points outside the 
global monopole's core can be expressed by 

dv 2 

ds 2 = rj^dx^dx" + — + r 2 (d6 2 + sin 2 6d<p 2 ) = g AB d A dx B , (1) 
or 

where rj^ = diag(—l , 1 , ... , 1) is the Minkowski metric, a 2 = 1 — k 2 ^, a parame- 
ter smaller than unity, and k related with the seven-dimensional Planck mass. In order 
to be more precise, in [13] the authors have obtained the solution to the Einstein equa- 
tions considering a general n— dimensional Minkowski brane worldsheet and a d > 3 
global monopole in the transverse extra dimensions. In this general case the metric is a 

2 2 

generalization of (CEJ) with fj, , v — 0,1... , n — 1 and a 2 = 1 — The solid angle 

associated with the 3— geometry above depends on the parameter a and reads Q = 4ir 2 a 2 , 
so smaller than the usual one. Consequently in this submanifold there is a solid angle 
deficit Ail = Att 2 k 2 t] 2 
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Composite topological defects have also been first studied in a four-dimensional space- 
time. Specifically a composite monopole, i.e., a system composed by a local and a global 
monopoles was analyzed in [TSJ UHl ED] (for composite strings see [21]). More recently, the 
composite monopole has been analysed in the braneworld scenario in |22j . 

In previous publication [23J, we have analyzed the vacuum polarization effects asso- 
ciated with a massless scalar quantum field in an (n + 3)— dimensional bulk spacetime 
which has the structure of a n— dimensional Minkowiski brane with a global monopole in 
the transverse three-dimensional sub-manifoldll] Specifically we calculated the renormal- 
ized vacuum expectation value of the square of the field, ($ 2 (x))R en , and have showed 
that this quantity depends crucially on the values attribute to n. We also analyzed the 
structure of the renormalized vacuum expectation value of the energy-momentum ten- 
sor, (Tab (x)) Ren.- In order to develop these investigations we calculated the respective 
Euclidean scalar Green function, G^\x,x'). 

Continuing in the same line of investigation, in this paper we shall analyze the po- 
larization effect of fermionic vacuum induced by a point-like three-dimensional global 
monopole embedded in a higher- dimensional bulk in a braneworld scenario. Specifically 
we shall consider that the bulk has its geometric structure given by the line element 
PI). Our main objective is to calculate the renormalized vacuum expectation value of the 
energy- momentum tensor, (Tab{%)) Ren.- Differently from the scalar case, here we obtain a 
simpler expression to the fermionic Green function and for the energy-momentum tensor, 
which is obtained in a closed form for an arbitrary parameter angle deficit. 

This paper is organized as follows: In section 2 we calculate the Euclidean fermionic 
propagator associated with a massless field in the background of a point-like global 
monopole transverse to a flat n— dimensional brane in a braneworld scenario, consid- 
ering n = 1, 2 and 3. In order to do that we write the general Dirac differential operator 
and the equation obeyed by the propagator. In section 3 we analyze this function in the 
coincidence limit and extract all divergencies from it in manifest form. In this way we pro- 
vide the explicit expressions to the components of the renormalized energy-momentum 
tensor, and analyze their behavior in some limiting cases. In section 4, we summarize 
our most important results. In this paper we use signature +2, and the definitions: 
Rp^s = <9 7 r^ 5 — R a /3 = Ray/3- We also use untis h—c — 1. 

2. Spinor Green Function 

Before to start the calculation of the spinor Green function in the six-dimensional global 
monopole spacetime, we shall first review briefly some important properties of the Dirac 
equation in a flat space. 

In a flat six-dimensional space, the Dirac matrices, T^ M \ are 8x8 matrices, which 
can be construct from the four-dimensional 4x4 ones [23] as shown below |25j: 

where 75 = Z7°7 1 7 2 7 3 , and / the 4x4 identity matrix. It can be easily verified that these 
matrices obey the Clifford algebra: {T^ M \ T^} = ~2^ M ^ N \ for M, N = 0, 1, ... 5. 

Although the physical interesting case corresponds to n = 4, we developed our formalism considering 
n as an arbitrary number 
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In these representation, the matrix is written as: 

,,.7! = .mOi.-ui ,...-,! _ / / 
...j - 1 Q _ ] . 



(3) 



This matrix has two chiral eigenstate defined as \&+ and Consequently any six- 

dimensional fermionic wave-function, can be decomposed in terms of its chiral compo- 
nents as^l 



Solutions for the Dirac equation, 



Mm 



(4) 



(5) 



with defined chirality can only be possible for M = 0. In this way for positive chirality, 
equation (jSJ) reduces to 

a^d iM) m + = , (6) 
being = (7^, —I) a set of 4 x 4 matrix; and for negative chirality it reduces to 

a iM) d {M )S>- = , (7) 

being now = (7^, 275, I). 

In order to write the Dirac equation in the six- dimensional global monopole spacetime, 
we shall choose the following coordinate system and basis tetrad: 



x A = (t, r,9, 0, x, y) = {x?, y a ) 
where \i = 0, 1, 2, 3 and a = 1, 2, and 



(8) 



e (M) 



( 1 

a sin 9 cos < 

a sin 9 sin < 
a cos 9 


y 



\ 

cos# cos0/r — sin0/rsin# 

cos#sin</>/r cos^/rsin^ 

-sin9/r 

10 

1 / 



(9) 



The Dirac equation for a massive field in the above coordinate system reads 

if^-M^ = 0, (10) 
with the covariant derivative operator given by 

^ = ef M) r (M) (^ + n A ), (11) 



2 It is worth to mention that the six-dimensional chirality does not correspond to four-dimensional 
chirality. Six-dimensional chiral wave-functions correspond to a four-dimensional wave- functions, which 
still contain two four-dimensional chiral components eigenstates of 7 s . 
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where IT^ is the spin connection, given in terms of the flat spacetime Dirac matrices by 



ii A — 4 1 e (M) e (N)C;A ■ U^J 







t -* A 

-(1 - a)E (8) ■ 












--(1 - a)sin6>£ (8) ■ 9 , 








-( 


, o 


I) Ml l) 



For the above basis tetrad, the only nonzero spin connections are: 



(13) 

where 

S(8) = {» " n *) , ( 14 ) 

being # and the standard unit vectors along the angular directions and a k the Pauli 
matrices. 

The fermionic propagator obeys the following differential equation: 

(if - M)S F (x,x') = — ^5 {6) (x -x')/ (8 ) , (15) 

where g = det^gAs) an d 1(8) is the 8x8 identity matrix. This propagator is a bispinor, 
i.e., it transforms as \l/ at point x and as \I/ at x' . 

If a bispinor D F (x',x) satisfies the differential equation below 



(d - M 2 - ~RJ D F (x,x') = --^=5^\x - x')/( 8 ) , 



(16) 



with generalized d'Alembertian operator given by 

□ = 9 MN V M V M = g MN {d M V N + U M V N - { S MN }V S ) , (17) 

the spinor Feynman propagator can be written as 

S F (x', x) = (if + M)D F (x', x) . (18) 

Now, after this brief review about the calculation of spinor Feynman propagator, let 
considering the six- dimensional global monopole spacetime ([Q), where the scalar curvature 
R — 2( ' 1 ~2 Q! - ■ Choosing the basis tetrad (Q, we obtain, after some intermediate steps, that 

1 .-.o o f ~n 2 _ \ Zj 2 _o _o (1 — (X 



!C = D- ^R = -d 2 + a 2 {d 2 . + -d r ) - — + dl + d 2 y - y -^-f^(l + ^ y L) , (19) 

being L the ordinary angular momentum operator. 

The system that we shall consider consists of a massless positive chiral fields. In this 
way Eq. ( [101) can be written in terms of a 4 x 4 matrix differential equation 

P® + = , (20) 
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with 

1 f -> -* \ n a. 



JP = 7°^ + a lr d r - - 7r (s • L + l) + 175^ - I0„ + ^7r , (21) 

being j r = f ■ j. 

The Feynman four-component propagator obeys the equation 

iJZ>S F (x, x') = -y=<5 (6) (x ~ , (22) 

and can be written in terms of the bispinor Gphy 

S F (x,x') =iflg F (x',x) , (23) 
where now Q F (x',x) obeys the 4x4 differential equation 

ICg F (x, x') = ^<5 (6) (x - x')I , (24) 



with 



« = -^ + ^(fif + ^)-^+«+^-2_^(i+f.r) . (25) 

The vacuum average value for the energy-momentum tensor can be expressed in terms 
of the Euclidean Green function. It is related with the ordinary Feynman Green function 
[26] by the relation Q F (r,r; r',r') = —iQ F (x,x'), where t = ir. In the following we shall 
consider the Euclidean Green function. 

In order to find a solution for the bispinor Q F (x,x'), we shall obtain the solution for 
the eigenvalue equation 

K E ®x{x) = -A 2 $ A (x) , (26) 

with A 2 > 0, so we can write 

>t, ,.M 

X 2 



WV)=E^. (27) 



Due to fact that our operator ( TT91) is self-adjoint, the set of its eigenfunctions constitutes 
a basis for the Hilbert space associated with four-component spinors. Moreover, because 
operator K, F is a parity even operator, its eigenfunctions present a defined parity, so the 
normalized eigenfunctions can be written as: 



tT ,->, ,., e ~ ikx fa ( JuAP^l^OA) \ , 



A 2 = k 2 + a 2 p\ (29) 



where kx = i] ab k a x b = k r + k x x + k y y, and J v represents the cylindrical Bessel function 
of order 

7 + 1/2 n 

v. = — - -rf , with n a = , a = 0, 1 . (30) 

a 2 
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These functions are specified by the set of quantum number (a, k a , p, j, rrij), where k a G 
(—00, 00), p G [0, 00), j = 1/2, 3/2, ... denotes the value of the total angular quantum 
number, rrij = —j, j determines its projection, a specifies two types of eigenfunctions 
with different parities corresponding to I = j — n a /2 being I the orbital quantum number. 
In fT28|) . ipj^j are ^ ne spinor spherical harmonics which are eigenfunctions of the operators 
L 2 and a ■ L as shown below: 

L\fX 3 = Kl + VtfX, (31) 
= + (32) 

with k(°> = — (I + 1) = — (j + 1/2) and = I = j + 1/2. Explicit form of above standard 
function are given in Ref. [23], for example. 

Although we have developed this formalism for a six-dimensional spacetime, it can be 
adapted to a four and five dimensional space. The reason resides in the representations for 
the Dirac matrices in these dimensions. For four dimensions a irreducible representation 
for the flat Dirac matrices is the well known 4x4 ones, so we may use a M = 7 M . 
As a consequence, in the corresponding analysis we have to discard the derivative with 
respect to the coordinates x and y in the operator /C. For a five dimensions a possible 
representation for the flat Dirac matrices is also 4x4 which can be given by o M = (7^, 175). 
In this case the coordinate y should be discarded. So on basis of these arguments it is 
possible to generalize the eigenfunction of the operator Ke as: 

~ (2vr)«/ 2 V r \in a J^(pr)r.<f<p<%(e,<l>) J ' 1 ) 

with n = 1, 2, 3. From now on, we shall use the above eigenfunction and make applica- 
tions for specific values attributed to n. 

Now we are in position to obtain the bispinor Qe, which is given by 

/co p poo 

ds I d n k / dp ^W^V) e_sA2 • (34) 

Finally substituting (PD into ([33]) with A 2 = k 2 + a 2 p 2 , we obtain with the help of [27J a 
closed expression to the Euclidean Green function: 

n-l 

G < r > r ) = 1 ( a \ 2 V ( F^oshu^) \ 

1 ' (2-Krr 1 )^ \ sinhtj F ijroj (coshu,n, 0') J ' 



(35) 



where 



n— 1 



F i>m .(cosh U ,C^) =Q„J^ , (36) 



being 



coshw = — — 7, > 1 , (Ax) 2 = f} ab (x - x'Yix - x'f , (37) 

2rr 
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and Cj^.(f), fi') = (pj^.{9,(j))(pj^.{9\(j)') a 2 x 2 matrix. is the associated Legendre 
function. We can express this function in terms of hypergeometric functions [27] by: 



Qt(coshu) = e^2^ ^ + + 3 J 2) J (1 _ e _ 2M)A+1/2 (sinh^ x 

Wa + 1/2, -A +1/2; Z/ + 3/2; jA^) ■ (38) 

In this analysis we identify the parameter v with v a — 1/2 and take A = (n — l)/2. So 
the relevant hypergeometric function is 

, ' n n — 2 

^ 1 77 , — ; v ° + 1 ; 



2 ' 2 ' ' u ' ' ' 1 - e 2u 

If n is a even number, this function becomes a polynomial of degree however being 
n an odd number, this function is an infinite series. 

Taking a = 1, u Q = j = 1 + 1/2 and V\ — j + 1 = Z + 1/2; however for the first 
index the angular quantum number I assumes crescent integer number starting from zero, 
i.e., I — 0,1, 2,..., while for the second index I — 1, 2, .... For I = the only 

term that contributes to (1351) in the summation is Q 2 (C)/ 2 1/2 + ^1/2 -1/2) ■ ^ or ^ — ^ 

both associated Legendre functions contribute, J2i>i Qi 2 m v^i,m + ^,m,)- Using the 
explicit expressions for the spinor spherical harmonics, and using the addition theorem 
for the spherical harmonics, it is possible to express fl35l) by: 

n-l 

fe^sH 1 _L(- * ) 2 ^(2/ + l)Q^(coahu)^(cos 7 )/. (39) 

Now taking n = 1, 2, 3 and using the sum of Legendre functions and polynomials 
we obtain standard expressions for the Green function. So: 

• For n — 1, we have 

Qe{x\x) = — 7Z7 ^I . (40) 

• For n = 2, we have 

-/ . (II) 

• For n = 3, we have 

<W*» = S5£^)i'- (42) 

After this application of our formalism, let us return to higher dimensional global 
monopole spacetime. In this space the fermionic Green function, Sp, can be given by 
applying the Dirac operator on the bispinor fl35l) according to (T23T) . 



1 


1 


4vr 2 (x' 


-x) 2 


1 


1 


8tt 2 (x 1 


— x) 3 


1 


1 



s 



3. Vacuum Average of the Energy- Momentum Ten- 
sor 

In this section we shall calculate in a explicit way, the renormalized vacuum expectation 
value (VEV) of the energy-momentum tensor, (T^Ren.- Because the metric tensor does 
not depend of any dimensional parameter, and also because we are working with natural 
units system, we can infer that the VEV of the energy-momentum tensor depends only 
on the radial coordinate r. Moreover, due to the nonvanishing of Riemann and Ricci 
tensor, and scalar curvature of the spacetime, this VEV also depends on the arbitrary 
mass scale, /x, introduced by the renormalization prescription. So, by dimensional analysis 
it is expected that 

(T£) Ren . = * (F§ + G A B Mjir/aj) , (43) 

where the tensors F B and G R depends only on the parameter a. Because the presence 
of the arbitrary cutoff scale, /i, there is an ambiguity in the definition of (|43p . Finally for 
a spacetime of odd dimension, i.e., for a even value of n, Gg = 0. Obviously the tensors 
are diagonal and due to the spherical symmetry of the problem we should have Fjj = F^ 

and G e e = G+. 

The renormalized VEV of the energy-momentum tensor must be conserved, 

V A {T*) Ren . = , (44) 
and provide the correct trace anomaly for spacetime of even dimension 



1 AT 

{TA)Ren - = J^f¥ Tr {a{n+3)/2) = (v /4¥r)"+3 ■ (45) 

As we shall see, taking into account these informations, it is possible to express all 
components Fg and G R in terms of the zero- zero ones, Fq and Gq, and the trace T 

Using the point-splitting procedure [26], the VEV of the energy-momentum tensor for 
this four-component spinor Feynman propagator, compatible with the eight- component 
one, has the following form: 

(T AB (x)) = ] lim Tr [a A (V B - V B ') + °b(V a - V A ')] S F (x, x') . (46) 

Because the dependence of the fermionic Green function on the time variable, the zero-zero 
component of the energy-momentum tensor reads: 

(Too(x)) = lim Tr lo d S F (x, x') , (47) 

x'—*x 

which can be expressed by 

(Too(x)) = -i lim d 2 t TrQ F {x\x) = - lim d 2 T TrQ E {x,x) . (48) 

x'^x x'^x 

In the obtainment of the above expression we have first taken in the Green function, 
fl35l) and fl36|) . the coincidence limit of the angular variable, Q = Q', and sum over rrij. 
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This makes this function proportional to the unit matrix I. So, only the term with time 
derivative in (1211) provides a nonvanishing contribution for the zero-zero component of the 
energy-momentum tensor. 

Now after these brief comments about general properties of the VEV of the energy- 
momentum tensor, we shall start the explicit calculation of this quantity, specifying the 
values for n. 



3..1 Case n = 1 

For n — 1, the bulk spacetime corresponds to the four-dimensional global monopole 
spacetime [29]. For this case the calculations of the vacuum polarization effects associated 
with massless two-component spinor field on this spacetime has been developed in |30j long 
time ago. More recently some Casimir densities associated with four-component massive 
fermionic field obeying the MIT bag boundary condition on, one spherical spherical shell, 
and two concentric spherical shells, in the global monopole spacetime have been analyzed 
in [31] and [32], respectively. 



3. .2 Case n 



The case n = 2 is a new one, it corresponds to a two-dimensional Minkowiski brane with 
a global monopole on the transverse sub-manifold. For this case the associated Legendre 
function in ( 1361) assumes a very simple expression 



Q 



1/2 



-1/2 



(coshu) 



(49) 



2 V sinh u 

Taking r = r' and Q — Q' into (I35p . summing over rrij and using the above expression 
for the Legendre function, it is possible to develop the sum over j, which is a geometric 
series, we obtain a closed formula for the Euclidean Green function: 

1 1 



Ge(x',x) 



a 



64ttV 3 sinh(u/2) sinh 2 (w/2a) ' 

with u = 2arcsinh(aAx/2r). 

The above Green function is divergent in the coincidence limit, Ax 
verify its singular behavior by expanding ( 1501 in power series of Ax: 



(50) 



0. We can 



Qe{x',x) 



a 



+ 



8vr 2 (Ax) 3 96vr 2 r 2 (Ax 
31a 6 - 21a 2 - 10 



+ 



a 



1920vr 2 r 4 



(Ax) 



483840vr 2 r 6 



(Ax) 3 + 0((Ax) 5 ) 



(51) 



In order to obtain a finite and well defined result to the VEV of the zero-zero compo- 
nent of the energy-momentum tensor, we must extract all divergent terms in the coinci- 
dence limit after applying the second time derivative. In order to do that in a manifest 
form, we subtract from the Green function the Hadamard one. In [2B] is given, for any 
dimensional spacetime, the general formal expression for the Hadamard function. For a 
five-dimensional spacetime it reads 

A 1/2 (x' x) 1 

2 V- — r [a {x,' , x) + ai(x', x)a(x', x) - a 2 (x', x)a 2 (x', x)\ .(52) 

107T V 2 O" [X , X) 



Gh{x , x) 
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For this manifold, A(x', x), the Van Vleck-Morette determinant, and the coefficients 
are given below: 



1 — a 2 1 — a 4 

A = 1 , a = I , ai = —5- I and a 2 = — ^Tri" 1 • ( 53 ) 

br z 60r 4 

The one-half of the geodesic distance for r' = r and ft' = Q, is a(x',x) = . So, we 
can see that the singular behavior for the Green function, Eq. (1511) . after taking its second 
time derivative has the same structure as given by the Hadamard function 



I -a 2 1 1 - a 4 



8vr 2 (Ax) 3 96n 2 r 2 (Ax) 19207r 2 r 4 



I . (54) 



So, on basis of this fact we can see that the renormalized VEV of zero-zero component 
of the energy-momentum vanishes, i.e., 

(T$(x)) Ren . = lim d 2 Tr [G E (x',x) - G H (x',x)] = . (55) 

x'— >x 

For an odd dimensional spacetime there is no trace anomaly, i.e., (T4 )# en . = 0. Writing 
Fg = (Fq, FjT, Fg, F%, F£), we have that the sum of all these components vanishes; 
moreover, the geometric structure of the brane section of this five- dimensional spacetime 
is Minkowiski-type; consequently the Green function and the Hadamard one depend on 
the variables on the two dimensional brane by Ax 2 = — (Ax°) 2 + (Aa; 4 ) 2 = (Ar) 2 + (Ax 4 ) 2 . 
By the definition of the VEV of the energy- momentum tensor, Eq. (I46I) . we have 



{Tt{x)) = lim Trd 2 x g E (x,x') . (56) 

x'— >x 

So, we can infer that (Tq(x)) R en , = (T£(x)) R en .. In this way we have F® = F£, being 
both zero. The conservation condition, V a(T^) R en . = 0, provides 3F r r = —Fg — Ft, and 

V A{Tg-) Ren, = 0, Fg = F%. So on basis of all these informations we conclude that all 
components of the tensor F^ are zero, consequently, 

(T^Ren. = • (57) 



3. .3 Case n = 3 

The case n = 3 corresponds to a flat three-dimensional brane transverse to a three- 
dimensional global monopole space. The respective Euclidean Green function is given in 
terms of an infinite sum of the associated Legendre function (5y tr _ 1 / 2 with C^.. Taking 
the angular coincide limit, Q = ft', the sum S below, given in fl35|) . can be expressed as: 

So = ^Qi _ 1/2 (coshn)C|°4(fi,fi) = + l)Qk_ 1 (co8h«)/ (2) (58) 

j,rrij l>0 

and for Si 

Si = J2^-y^oshu)C^.(n,n) = i-^/Q| (cosh u)/ {2) . (59) 

j,mj l>l 



ll 



Using Ql(z) = {z 2 — i)V2 anc j integral representation below for the Legendre 
function Q u 

-{v+l/2)t 



Q v (cosh.u) 



dt 



\/2 Ju V cosh t — cosh u 
it is possible to develop the sums over I in (T3ST) and floTJj) getting: 



and 



So (cosh u) 



5*i (cosh u) 



1(2) sinh u d 
dz 



1(2) sinh u d 
1QtiV2 dz 



dt e 1 / 2 



arccoshz Vcosh t - z sinh 2 (t/2a) 
dt e"*/ 2 1 



z=cosh u 



z=cosh u 



(60) 



(61) 



(62) 



by 



/arccoshz Vcosh t - z sinh 2 (£/2a) 
Substituting the sum Sq + Si into ( |35|) . we can express the Euclidean Green function 



Qe{x',x) 



a 2 d 



32?r 3 (r'r) 2 dz 



dy 



— VV + 1 - z sinh 2 



2 / arcsin 



2=COsh M 



(63) 



where we have introduced a new variable t := 2 arcsinh(y/\/2). 

The VEV of the zero-zero component of the energy-momentum tensor can be formally 
given by substituting (I6"3"|) into (14"8"I) . Because this procedure provides a divergent result, 
we must renormalize it by extracting all the divergent terms. In what follows we shall 
develop a procedure to extract the divercencies in a manifest form. Let us consider the 
integral inside the bracket of (163]) and write it as: 



UZ) 



+ 



h(z)+I 2 (z) 
dy 



dy 



~ VV + 1 - z sinh 2 



2 / arcsin 



1 a/?/ 2 + 1-2 s i n h 2 ^ arcsinhCy/y'S) 



(64) 



Subtracting and adding into the integrand of I\ the first four terms of power series of the 
function 



sinh 

we get 



2 / arcsinh(y/v / 2) A y^ 
« J 



2a 2 (a 2 - V 
+ - — - — - 



fa 4 - 1^ 



30a 2 



> , (31^-210^-10) 

3780a 4 y V ; 



^(^ = i[ m ( z )+ir 9 {z) 



(66) 



with 

r fin I \ 

1[ (z) 



dy 



=1 y/y 2 + 1 - z 

(31a 6 - 21a 2 - 10) 
3780a 4 



2a 2 



fa 2 - 1) (a 4 - 1) 



sin h 2 ( arcsinh(y/y^) \ y2 



+ 



30a 2 



-y 



(67) 
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and 

Il in \z) 



dy 



2a 2 {a 2 - 1) (a 4 - 
~y*~ + 3 30a 2 



1) 2 (31a 6 -21a 2 



3780a 4 



.(61 



As we shall see, l( m together with J 2 , provide a finite contribution to the VEV. All 
the divergences are contained into I{ m9 . 

Taking the derivative of 7* m9 with respect to z and substituting z = coshu = 1 + 

with a = ^y~, we expand the result in power of which survives after taking the second 
derivative and the coincidence limit. Doing this procedure we have: 



dz 



2r 4 



a 2 a 2 



a x r 2 a 2 

— 7; In 

2a 2 a 4 



era 



{-Aa 2 - 10a 3 + 6a 2 - 3a x - 6a 3 In fe)) ™ 2 



16r 2 



+ 



with 



fa 2 - r 



a 



«2 



1) ,_ (31a 6 -21a 2 -10) 
and a 3 = — 



3780a 4 



(69) 



(70) 



3 ' 30a 2 

being last three coefficients of the expansion fl65|) . 

We can see that the contribution to the Euclidean Green function given by I^™ 9 has 
the same structure as the Hadamard function for a six- dimensional spacetime, 



Gh{x', x) 



A 1 / 2 

16/T 3 



a tzi 1 



2(7 



) In 



~4~ 



(71) 



where /1 is an arbitrary energy scale [28]. A, do, Qi and a 2 are same as given in (|53|) for 
this six-dimensional global monopole spacetime. 

In [331 [34] , an explicit expression for the coefficient a 3 is provided for a general second 
order differential operator D 2 +X, Dm being the covariant derivative including gauge field 
and X an arbitrary scalar function. For this six- dimensional global monopole spacetime, 
by using the computer program GRTensorll, we found, considering X = —jRI, 

31a 6 - 21a 2 
a 3 = - 



10 



-I . 



2520r 6 

The renormalized bispinor given by 

GRen.(x',x) = Ge{x',x) - G H (x',x) 

with 



Ge(x'x) 



a 



32vr 3 r 4 



sing 



dz 



a) + dl{ 



fin. 



dz 



er) dl 2 (<j) 

dz 



(72) 
(73) 

(74) 



and Gh(x',x) given by (ITTj) . provides 



QRen.{x',x) = I 
- I 



a 



647r 3 r 4 
2 



a 2 



3a 3 a 2 cx 
2r 2 



ln(^ 

a 



[—4a 2 — 10a 3 + 6a 2 — 3ai]aV 
5127r 3 r 6 



a 



32vr 3 r 4 



d/f n ((7) d/ 2 ((7) 



(75) 
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Now we are in position to obtain the renormalized VEV of the zero-zero component 
of the energy momentum tensor. It is given by: 

(T$(x)) Ren . = lim d 2 T Tr g Ren {x\x) . (76) 

Wald in |35j, has proved that in order to obtain a energy-momentum tensor which 
obeys the conservation condition law (1441) and correct trace anomaly (T^oT) . an additional 
contribution must be considered. In a six- dimensional spacetime, this terms reads: 

<^Tro 3 • (77) 



384vr 3 



So on basis on this fact we explicitly present in Appendix |A] the steps needed to obtain 
Fq and Gq. Particularment with respect to Fq we provide an integral expression for it, 
and for Gq a closed expression. Briefly speaking, they are given by substituting ( |75l) into 
( 1761) and taking into account (J77J. These components are: 

F = _ 6a 4 / OS /i(x) _ fo 4 / ^ h{x) + g a 4- 3 (7g) 
.In X ./i X 



with 



fix) = 1 2a2 (o: 2 — 1) (a 4 — 1) ^2 

sinh 2 (^^kl^l) x 2 3 30a' 



\ a I 

(31a 6 - 21a 2 - 10) , 



3780a 4 



(79) 



. . 1 2a 2 (a 2 - 1) (a 4 - 1) , 



sinh 2 



and 



Gq = -r b Tr a 3 



630 



In (1751) we have used the fact that (r 6 /6)Tra 3 = a 4 a 3 . 

Unfortunately it is not possible to provide analitycal results to the integrals above 
for a general value of a. Their dependence on the parameter a can only be provided 
numerically. Our numerical results to these integrals are exhibited in figure [TJ Also 
it is possible to provide the approximate behavior for the integrals, l\ = / jir/i and 
h = J] 00 tt/2> f° r specific limit of this parameter: 

• For large solid angle deficit (a << 1), 

/ x 1 ln(a) . , 

Ha) « --L,. (83) 
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For small solid angle deficit (|a — 1| << 1), 

h(a) « -0.0052(a - 1) 
J 2 (a) « 0.0275(a- 1) . 

For large angle solid excess (a >> 1), 

7i(a) w -0.0025a 2 
a 2 

V ; 60 



54) 
55) 



(86) 
(87) 




0.04- 






0.02- 


^ — alpha 

0.5 IS 


2.5 




4.02- 






4.04- 


/ 




-0.06- 







Figure 1: These graphs represent the dependence of the two integrals in (1751) as a function 
of the parameter a in a specific interval. The left panel corresponds to integral associated 
with the function fi(x)/x 5 , and the right panel corresponds to integral associated with 
the function f 2 (x)/x b . 

Having found Fq and Gq the others components of the renormalized VEV of the 
energy-momentum tensor can be expressed in terms of them. Due to the Minkowiski 
structure of the brane section of this spacetime, the Green function and also the Hadamard 
one depend on the variables on the brane by (Ax) 2 
hand, by using (1461) we can verify that 



Vab{X 



x') a (x 



X 



t\b. 



on the other 



(T a a (x)) = lim Tr d 2 a Q E (x',x) , for a = 4, 5 . (88) 

So we may conclude that the renormalized VEV of these components are equal to zero-zero 
one. Consequently F$ = 



Fl and G° 



G\ = G\. Besides the trace anomaly 



(Ti(x)) 



Ren. 



64vr 3 



Traa 



T 



167r 3 r 6 



with 



T 



31a 6 - 21a 2 
2520 



10 



(89) 



(90) 
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give us F£ = 4T and = 0. The conservation conditions, V a(T^) R en . = 0, and 
V A (T r A ) Ren . = 0, provide, respectively, T e e = T$ and 4F r r - G r r + 2F e B = 0. So after a 
simple calculation we can express: 

Fb = diag(F °, F ° + G°/3-4T/3, 8T/3 - 2F ° - G°/6, 

8T/3 - 2F ° - G°/6, F °, F °) (91) 

and 

G£ = G°diag(l, 1, -2, -2, 1, 1) . (92) 

4. Concluding Remarks 

In this paper we have considered the vacuum polarization effect associated with fermionic 
field induced by a global monopole in the braneworld context. By this scenario our 
Universe is described by a (n — 1)— flat brane transverse to a point-like three-dimensional 
global monopole manifold. Two specific spacetime have been explicitly considered, the 
cases with n = 2 and 3. In fact our first motivation was to consider a six-dimensional 
bulk. However, because we have considered a defined positive chiral field, the fermionic 
field can be described by a four-component representation; so, in this way, we could 
extend the formalism to include a five- dimensional bulk. Due to the fact that the global 
monopole lives in a three-dimensional manifold, it was possible to express the fermionic 
eigenf unctions in terms of spinor spherical harmonics, v 9 ^! • Differently from the scalar 
case analyzed in [23], the effective total angular momentum have simple expressions: 
u = (Z + I) /a - 1/2 and v x = l/a + l/2@ 

The main objective of this paper was to obtain the renormalized vacuum expectation 
value of the energy-momentum tensor, (T^ )■ In order to do that we have explicitly con- 
structed the fermionic Green function, which was expressed in terms of a linear differential 
operator acting on a bispinor. Because all components of this tensor can be related by 
the conservation condition and the correct trace anomaly, we needed only to calculate 
the zero-zero component of this tensor. By our results we could verify that for a five- 
dimensional bulk (T4 )_R en . = 0. @ However, for a six-dimensional bulk, a non- vanishing 
result have been obtained. This result depends on the inverse of the sixth order power of 
the distance from the point to the monopole's core considered on the brane, and also on 
the arbitrary energy scale fi: 

(T °(x)) Ren . = (F V) + G° (a) Hfir/a)) . (93) 

The expressions found for F® involve a long calculation and some details related with 
this calculation is presented in Appendix |A3 By our result this component is given in 
terms of a closed term and two more integral expressions. For these integrals we provided 
numerically, in Fig. [U their behavior as function of a. Fortunately we were able to 

3 For the scalar case the effective orbital angular quantum number assumes a simple form only for a 
curvature coupling parameter £ = 1/8. 

4 A vanishing result for the renormalized vacuum expectation value of the square of the field, ($ 2 ), 
have also been obtained for a scalar field in this five-dimensional bulk in a specific approximated result 



16 



provide a closed expression for Gq, which by its turn depends on the coefficient 03 of the 
Hadamard function. The general expression for this coefficient, for the fermionic case, 
presents 43 terms, so its final result requires a long and careful calculation, even using 
the specific computer program GRTensorll. Because fl93|) is given in terms of two distinct 
contributions, there exist a radius f which (T® (x))n en . vanishes. In figure [2] we present 
the dependence of the dimensionless distance [if jot as a function of a. 




Figure 2: This graph represents the dependence of the physical distance fir/ a, that 
(Tq (x))n en . vanishes, with a. 

Finally we would like to mention that the model analysed here presents monopole as 
a point-like object having its core on the flat brane; so the influence due to the fermionic 
quantum field on the brane can be evaluated in the region near the monopole's core. 
However, by this model the renormalized VEV of the energy-momentum tensor is di- 
vergent at the monopole's core. The problem of the singularity of vacuum polarization 
effects involving a global monopole, and topological deffects in general, can be avoided by 
considering a more realistic model to the monopole, i.e., considering a inner structure to 
its core. A simplified model for the monopole core has been first presented in [36]. The 
vacuum polarization effects due to a massless scalar field in the region outside this model 
have been investigated in [27]. More recently the analysis of vacuum polarization effect 
associated with quantum bosonic and fermionic fields, in the global monopole spacetime 
with a general spherically symmetric inner structure, have been developed in |38j and 
|39j . respectively. In these analysis the asymptotic behavior of the core- induced vacuum 
densities are investigated at large and small distance from the core, and for small a large 
angle solid deficit. 
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A Explicit Calculation of Fq and G® 

In order to find the complete expression for (!75|) . it is necessary to obtain the derivative 
of l( m and I2 with respect to z. Our first steps in this direction is to introduced a new 
variable b = \J z — 1. So we may write: 



T fin/ \ 

I{ {z) 



dx 



with 



Vx 2 - b 2 



2a 2 



(94) 



a 



sinh 



2 / arcsin 



h(x/V2) \ 



+ 



(a 4 



30a 2 



-x 



(31a 6 - 21a 2 - 10) 
3780a 4 



-x 



Consequently 



dx 



lj in {a) _ ]_d_ f 1 

dz ~ 2b db J b Vx 2 - b 2 



fi(x) ■ 



(95) 



(96) 



Because the integrand of l{ m is divergent at the point x = b, we have to change the 
variable x — > bx before applying the Leibnitz formula for derivative of an integral. So, 
this derivative can be written as 



dx 



dz 2 J b ^x 2 -b 2 2VT^¥ 
where g(x) = and b = The derivative of I2 is simple and reads 



(97) 



dz 



dx 



2 J 1 (X 2 - b 2 ) 3 / 2 ginh 2 ^ arcsinh(x/V2) ^ 



(9c 



Using these results we can express the renormalized Green function, Eq. ([731) . as: 



IReri. 



(x , x) 



a 



647r 3 r 4 
,2 



a 2 



a 



64vr 3 r 4 

„2 



- / 



a 



647r 3 r 



3a3a 2 cr 
" 2r 2 

da; 

6 \Jx 2 — b 2 
dx 



ln[^ 



a 



-4a 2 — 10a3 + 602 — 3ai]a 4 cx 
512vr 3 r 6 



Vl^J 
1 



52)3/2 



sinh 



2 / arcsin 



inh^x/\/2A 



(99) 



The next step, according to (176]) . is to obtain the second derivative with respect to 
the Euclidean time. Taking the coincidence limits x' = x and y' = y, using the fact that 
d 2 = (a 2 /2r 2 )d 2 , and adopting the same procedure as explained above we can write 



db 2 



dx 



y/x 2 - b 2 



9{x) 



dx 



, Vx 2 - b 2 
g(l)b 2 



„ ( \ 9'{x) g(x) 
g (x) + — ^ 



x 



52)3/2 



(100) 
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Taking the limit b — > we obtain 



' =flr(a:) -> / dx . (101) 



db 2 J b 7 7 ^ 3 

We can see that the integral above is finite because for small values of x, g(x) — > 
_ (289a -i68a -^looa ~ 2l - > x 4 + 0(x 6 ). In the limit 6 — > 0, the other relevant terms in the 
renormalized Green function above provide 

1 3 (102) 

1 . (103) 



db 2 (x — 6 2 ) 3 / 2 x 5 
d 2 1 



db 2 y^W 

So taking into account all these results and including the additional term (ITT)) we have: 
a 6 [6a2 — 3ai — 10a 3 ]a 4 3a 4 <ir 



327T 3 r 6 1287T 3 r 6 327T 3 r 6 ./l £ 5 sinh 2 f arcsinh(x/V2) 



3a 4 f 1 dx j, , s i 3a3« 4 , ( \it\ i , 



+ 32^ /„ + 3W ln - 64^" " 3 ' (1 ° 4) 



which can be written as shown below: 



3a 4 f 1 dx 3a 4 /"°° dx 

<T 00 (*)w = ^,y o ^/iw + ^l ^) 

3a 4 _ 3a 3 a 4 //ir\ 

i a3 + T^-r^ln — • ( 105 ) 



32vr 3 r 6 32vr 3 r 6 V a 
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